From a single-band metal to a high-temperature 
superconductor via two thermal phase transitions 
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The nature of the pseudogap phase of cuprate high-temperature superconductors 
is a major unsolved problem in condensed matter physics. We studied the com- 
mencement of the pseudogap state at temperature T* using three different tech- 
niques (angle-resolved photoemission spectroscopy, polar Kerr effect, and time- 
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resolved reflectivity) on the same optimally-doped Bi2201 crystals. We observed 
the coincident, abrupt onset at T* of a particle-hole asymmetric antinodal gap in 
the electronic spectrum, a Kerr rotation in the reflected light polarization, and a 
change in the ultrafast relaxational dynamics, consistent with a phase transition. 
Upon further cooling, spectroscopic signatures of superconductivity begin to grow 
close to the superconducting transition temperature (T c ), entangled in an energy- 
momentum dependent fashion with the pre-existing pseudogap features, ushering 
in a ground state with coexisting orders. 

As complex oxides, cuprate superconductors belong to a class of materials which exhibit many 
broken-symmetry states; unravelling the relationship between superconductivity in the cuprates and 
other possible broken-symmetry states has been a major challenge of condensed matter physics. A 
possibly related issue concerns the nature of the pseudogap in the cuprates and its relationship with 
superconductivity. Angle-resolved photoemission spectroscopy (ARPES) studies have shown that 
the pseudogap develops below a temperature T* near the Brillouin zone boundary while preserving 
a gapless Fermi arc near the zone diagonal (/). A key issue is the extent to which the pseudogap 
is a consequence of superconducting fluctuations (2-5), which should exhibit a rough particle-hole 
symmetry, or another form of (incipient) order (6-12), which typically should induce particle-hole 
asymmetric spectral changes. Candidate orders include various forms of density wave, nematic or 
unconventional magnetic orders that break different combinations of lattice translational (6-8, 13- 
19), rotational (6, 9, 15, 17, 20-22), and time-reversal (7, 9, 23-26) symmetries. 

We have focused on crystals of nearly optimally-doped (OP) Pb .55Bii. 5 Sri. 6 Lao.4Cu06+5 (Pb- 
Bi2201, T c = 38 K, T* = 132 ± 8 K) (27), and combined the ARPES measurements of the 
evolution of the band structure over a wide range of temperature, momentum and energy, with 
high-precision measurements of the polar Kerr effect (PKE) and time-resolved reflectivity (TRR). 
Bi2201 was chosen to avoid the complications resulting from bilayer splitting and strong antinodal 
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bosonic mode coupling inherent to Bi 2 Sr 2 CaCu20s+5 (Bi2212) (1). Whereas ARPES is a surface 
probe, PKE enables us to monitor a bulk, thermodynamic (via the fluctuation-dissipation theorem) 
property which has proven (28) to be a sensitive probe of the onset of a broken-symmetry state, and 
TRR gives complementary information on the bulk, near-equilibrium dynamics of the system. 

We will first analyze our ARPES data collected in different temperature regions. Above T*, Pb- 
Bi2201 has a simple one-band band structure (right side of Fig. [[]). For each cut in momentum space 
perpendicular to T-M [(0, 0)-(7r, 0)] (C1-C7 in Fig. [[]), the only distinct feature in the corresponding 
Fermi-function-divided (27) energy distribution curves (EDCs) is a maximum (red circles in Fig. |2l 
A-G). As a function of the y component of the wave vector (k y ), the maxima have an approximately 
parabolic dispersion for each cut (red circles in Fig. |2] O-U); the band bottom lies on the T-M axis, 
and the dispersion crosses the Fermi level (E F ) at two momenta, k F (k F1 and k F2 ). The band- 
bottom binding energy monotonically decreases from near T to M (Fig. [2l O-U). We take the 
Fermi-level crossings of this single band to define the Fermi surface. Despite the simplicity of the 
electronic structure above T*, the width and energy-dependent broadening of the EDC maximum 
features, along with the familiar strange metal behavior seen in transport, imply that the system is 
not well described as a Fermi liquid. 

We now turn to the temperature region below T c . Here, the entire Fermi surface is gapped except 
at the nodal points (k F lying on the zone diagonal). In the nodal region, consistent with previous 
reports (4, 5, 11, 12), a ci-wave-like gap along the Fermi surface is observed which we quantify as 
the energy position of the EDC maximum (blue circles) at k F (Fig. [2J L-N). This maximum is still 
the only identifiable feature in the EDC. By comparing the EDCs in Fig. [2l E-G with those in Fig. 
[H L-N, we see that the peaks of EDCs near k F are much sharper below T c than above T*, however 
(perhaps surprisingly) the peaks well away from k F appear broader but with larger experimental 
uncertainties (also see Fig. |2] V-W). 

Away from the nodal region, the dispersion along each cut rises to a minimum binding energy 
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and then bends back (Fig. [2l H-K). These back-bendings (black arrows in Fig. |2] O-S) occur at 
momenta k G (kai and ka2), which are increasingly separated from the Fermi surface (compare blue 
and red squares on the left side of Fig. [[) towards the antinodes (k F lying on the zone boundary). 
Note that, for a superconducting gap, as a consequence of the particle-hole symmetry, one would 
expect kc — kp (Fig. S6), as is the case in the nodal region. The substantial kc — kp misalignment, 
previously found in a single antinodal cut perpendicular to T-M, has been interpreted as a signa- 
ture of non-particle-hole symmetric order (29). Our observation here complements that finding by 
providing an entire momentum-space picture. 

The width and shape of the EDC near k F change fairly abruptly as a function of position along 
the Fermi surface, in contrast with the smooth evolution seen above T* (compare magenta curves 
in Fig. |2l H-N and A-G). Strikingly, the antinodal EDC maxima below T c are broader than those 
above T* (Figs. |2l V-W and S7I). Such broadening of the antinodal spectra with decreasing tem- 
perature, as reported previously (29), is likely intrinsic, as the expected sharpening is observed 
simultaneously in the nodal spectra at low energy. 

Around the M point, both the EDC line shape and the dispersion of the EDC maxima are more 
complex than their counterparts above T*. The dispersion of the EDC maxima has two separate 
branches, one at relatively low energy which shows back-bendings at k G and the other at higher 
energy around the band bottom on the T-M axis. As one moves from cuts CI to C7, the apparent 
discontinuity between the two branches near the antinode vanishes and the two branches merge 
(Fig. |2l O-U). Both k G (blue squares in Fig. [T) and the band-bottom energy (blue circles in Fig. 
|2ftV) evolve smoothly, even as the dispersion at intermediate energies changes dramatically. 

Along with the changes in the dispersion of the EDC maxima below T c , a well-defined shoulder 
feature (green dots in Fig. |2] H-K and W) (27) emerges at low energy in the EDCs. This feature 
exhibits little dispersion either along T-M (Fig. [2]W) or perpendicular to it (Fig. |2] H-J). Unlike the 
EDC maxima, its dispersion does not continue towards the zone center (Fig. |2M0 but instead loses 
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its definition away from the vicinity of the M point (magenta- shaded region in Fig. [T). 

The changes in the spectral function on cooling do not occur smoothly; instead, abrupt changes 
in the thermal evolution occur in the neighborhoods of T* and T c . 

A detailed temperature-dependent study of an antinodal cut (CI in Fig. [D shows that the dis- 
persion of the EDC maxima (blue circles in Fig. SI, A-E) exhibit a transformation which begins 
at a temperature equal (within experimental uncertainty) to the reported values of T* in the lit- 
erature (27). Specifically, the energies of the EDC maxima at the band bottom and at kp are 
temperature independent above T* but begin shifting to higher binding energies below T* (Fig. 
S1F). The energy position at k F is a measure of the pseudogap, and can be defined as a spectral 
order parameter that becomes non-zero below T*. Moreover, at temperatures below T* but well 
above T c , one can already see the k G — k F misalignment developing (right inset of Fig. |3]). These 
observations are consistent with our previous report on a similar cut in an OP Pb-Bi2201 sample 
prepared under a different post- annealing condition (29). 

It would be natural (as discussed below) to associate the abrupt change in nature of the ARPES 
spectra with a transition to a broken- symmetry state below T*. To test this idea, we performed 
PKE measurements on the same crystals with finely-spaced temperature steps (27). The results 
(Fig. [3]) are clearly suggestive of a slightly -rounded phase transition at T* below which a finite 
Kerr rotation emerges. A similar transition in the Kerr rotation was previously observed in high- 
quality YBa 2 Cu 3 6+;r (YBCO) crystals and was suggested to be linked to a broken- symmetry state 
that is not necessarily magnetic in nature (25). The smallness of the previously observed Kerr 
rotation in YBCO suggested that it might not reflect the primary order. The present data again 
show a small Kerr rotation; however, the strong correspondence between the ARPES and PKE data 
in Fig. |3l allows us to conjecture a phase transition at T*, and also to corroborate the interpretation 
of previously published data on YBCO in terms of a similar transition. 

We also studied TRR on the same Pb-Bi2201 crystals (27). At temperatures between T c and T*, 
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we observed a negative signal with a time-resolution-limited turn-on followed by decay on the order 
of 100 femtoseconds (Fig. S3). Additionally, we observed a positive signal that emerges below T c 
and decays on a picosecond timescale. These observations are consistent with previous studies of 
various cuprate families that revealed different dynamics above and below T c (30), and associated 
the former with the pseudogap state (31-36). The striking aspects of the TRR results reported here 
are the clarity of the onset at T* and the direct correspondence with ARPES and PKE. As shown in 
the left inset of Fig. [3l the magnitude of the negative TRR signal tracks the ARPES and PKE data 
quite well. This added correspondence further supports the existence of a phase transition at T* to 
the pseudogap state with near-equilibrium dynamics distinct from superconductivity. 

Upon cooling below T c , it is observed in ARPES that the shoulder on the low-energy side of 
the EDC maximum appears to develop somewhat above T c but well below T*, and it grows truly 
distinct only below T c (compare Figs. 51 A-B and SI, A-E). In order to obtain a clearer view of 
the structure of this feature, we show, in the insets of Fig. |4l E-G, the spectra at low temperatures 
divided by the one at 60K - a procedure that converts the shoulder into a small peak. In Fig. 
we adopt the same procedure, but dividing by the 40K spectra instead, and in Fig. S2 we have 
subtracted off an approximate background from the EDCs. All three methods of analysis produce 
qualitatively similar results, with a small peak at a position that does not disperse appreciably along 
the exemplary antinodal Cut CI. In all cases, the energy and width of the peak do not change sig- 
nificantly with increasing temperature, while the peak intensity is strongly temperature-dependent; 
the peak becomes undetectable at a temperature which, while greater than T c , is nowhere near T*. 
This behavior of the processed data is highly reminiscent of the behavior of the superconducting 
coherence peak in Bi2212 (37), although that peak is visible in the raw data (Fig. |4j3). 

The strong and analogous temperature dependences of the ARPES, PKE and TRR data seen 
between T c and T* are most naturally understood if T* is associated with a phase transition to a non- 
superconducting broken-symmetry state (27). While superconducting fluctuations are observable 
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above T c (Figs. HI E-G and S2E), they appear to peter out at too low temperatures (being negligible 
already at 60 K) to play a central role in this higher temperature transition. To complete this picture, 
it is necessary to identify the nature of the broken-symmetry state, and to relate it to the apparently 
similar electronic changes that occur below T* in other cuprates (27). 

Below T c , the nodal arc is gapped with a d-wave-like structure suggestive of a dominantly 
superconducting origin (38). In contrast, in the antinodal region, rather than one order being domi- 
nant, or the two gaps of both orders adding in quadrature, the spectral function develops a complex 
structure with two energy scales below Ep of mixed origin, a larger one being primarily associated 
with the pseudogap order, and a smaller one with the superconducting order. 

To see what can be learned about the pseudogap order from ARPES, we have used a simple 
mean-field model (27) to compute the expected changes to the band structure induced by various 
forms of density wave (6-8, 13-19) (Fig. S7) or nematic order (6, 15, 17, 20-22) (Fig. S8) coex- 
isting with c/-wave superconductivity. We note that some (but certainly not all) key aspects of our 
experimental observations can be qualitatively reproduced, regardless of which of these orders is 
assumed (27). An essential feature of both the experiment and the fits is the comparable sizes of the 
superconducting gap and the pseudogap. This implies that the two orders may have a more intimate 
connection than just competing orders, such as seen in 2i/-NbSe 2 where the charge density wave 
gap is at least 3 times the superconducting gap (40). 
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Figure 1: Fermi surface maps measured below T c at 10 K (left) and above T* at 172 K (right) in 
the same momentum- space region (flipped for display). Dashed white lines labeled C1-C7 depict 
the cuts along which the EDCs shown in Fig. |2l A-N were measured. Purple dots labeled P1-P16 
along M-r indicate momenta where EDCs in Fig. [21 V-W) were measured. Red and blue squares 
on the left indicate momenta of the Fermi crossing k F (kpi and k F2 in Fig. [2l A-G) at 172 K and 
back-bending k G (black arrows in Fig. |2] O-S) at 10 K of the dispersion of the EDC maximum 
along cuts C1-C7. Red and blue circles on the right indicate momenta of identifiable peaks in the 
momentum distribution curves (measured along cuts parallel to Cut C7) at Ep at 172 K and 10 K, 
respectively. The solid red curves are a guide to the eye for the red squares and circles, whereas 
the dashed blue curve is the guide for the blue squares, together they show an increased k G — k F 
misalignment going away from the nodal towards the antinodal region. The magenta- shaded region 
is approximately where multiple EDC features are found at 10 K (Figs. [2]W and S2F). 
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Figure 2: (A)-(G) and (H)-(N) Selected EDCs at 172 K and at 10 K, respectively, for cuts C1-C7 
approximately perpendicular to T-M (Fig. [T]). Each EDC corresponds to a white point in the cuts 
in Fig. \T\ EDCs in magenta and orange are located close to k F . (O)-(U) Dispersions of the EDC 
features in (A)-(N) for cuts C1-C7. For each dispersion curve, every other symbol corresponds to 
an EDC in (A)-(N). Error bars are estimated based on the sharpness of features, to be ±3 meV 
minimum and ±8 meV maximum [examples shown in (O)] based on various EDC analyses (27). 
(V)-(W) EDCs at momenta PI -PI 6 along M-T (Fig. [B at 172 K and 10 K, respectively. Circles 
denote the EDC shoulder feature (solid green), the EDC maximum feature at 10 K (blue) and at 
172 K (red). 
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Figure 3: Temperature dependence of Kerr rotation (9k) measured by PKE, in comparison with that 
of the binding energy position of the EDC maximum at kp given by ARPES (reproduced from Fig. 
S1F and Ref. (29)). ARPES results are normalized to the 80 K values (free from the interference 
of fluctuating superconductivity). The dashed black curve is a guide to the eye for the PKE data, 
showing a mean- field-like critical behavior close to T* (see additional discussion (27)). Left inset: 
Temperature dependence of the transient reflectivity change measured by TRR (right axis). The 
dashed black curve (left axis) is reproduced from the main panel. Error bars (if invisible) are smaller 
than the symbol size. Right inset: The dispersion of the EDC maximum at various temperatures 
above T c , summarizing the results of Figs. UK, SI, A-E and|4]A All data were taken on samples 
from the same growth and annealing batch, except those reproduced from Ref. (29) on differently 
annealed samples. 
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Figure 4: (A) and (B) Selected EDCs at 40 K and 22 K along Cut CI (Fig. See Fig. A and 
H for data at 172 K and 10 K, Fig. SI, A-E for other intermediate temperatures. (C) Antinodal 
EDCs at 10 K after dividing by the 40 K counterparts, covering the momentum range indicated 
by the grey bar in (H), in comparison with those in (D) taken in a similar range at 30 K on an OP 
Bi2212 sample. Non-dispersive peaks are seen in both cases despite different sharpness and energy 
positions. (E)-(G) EDCs at different fixed momenta [specified in (A) and (H)] and temperatures 
around T c . The counterintuitive increase of the antinodal gap, defined by the energy position of the 
EDC maximum in (F) and (G), with temperature rising above T c cannot be understood with a single 
energy scale assumed. Insets: Corresponding EDCs divided by the 60 K counterpart, showing the 
peaks losing definition above T c (Fig. S2E). (H) Summary for the dispersions of related EDC 
features across and below T c . Vertical arrows specify momenta M, k,F2 at 172 K and ko2 at 10 K. 
Apparent asymmetry of the dispersions across M is due to a finite deviation of the cut from the 
high- symmetry direction and a subtle balance of spectral weight between different features in the 
EDC. All EDC features and error bars are similarly determined as in Fig. [2] 
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